Angular momentum exchange between coherent hght and matter fields 
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Full, three dimensional, time-dependent simulations are presented demonstrating the quantized 
transfer of angular momentum to a Bose-Einstein condensate from a laser carrying orbital angular 
momentum in a Laguerre- Gaussian mode. The process is described in terms of coherent Bragg 
scattering of atoms from a chiral optical lattice. The transfer efficiency and the angular momentum 
content of the output coupled vortex state are analyzed and compared with a recent experiment. 
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Photons possess both translational and spin-angular 
momentum, where the latter is associated with differ- 
ent polarization states of the light wave. Light may also 
be prepared in states with well-defined orbital angular 
momentum as demonstrated through the generation of 
Laguerre- Gaussian (LG) beams [JJ l2j. Coherent con- 
trol of the interaction between light and matter facili- 
tates macroscopic transfer of those quantities from one 
medium to another, opening up exciting and novel possi- 
bilities for quantum information engineering such as stor- 
ing and transferring light in ensembles of atoms [3l |4] . 

In a recent experiment, orbital angular momentum 
was controllably transferred in discrete quanta from laser 
photons to a Bose-Einstein condensate (BEG) of atoms 
[5] . In this type of experiment the creation of a quantized 
vortex via coherent transfer of angular momentum from 
photons to atoms is direct and therefore does not rely 
on mechanical rotation or subsequent (thermo) dynamical 
equilibration. Furthermore, this process provides a clean 
way of preparing topologically-distinct quantum states 
and enables the creation of genuinely multiply-charged 
vortices of any integer winding number. It also serves as 
a prototype for coherent generation of more complicated 
states by employing the spatial phase of the wavefunc- 
tion. Further development may produce a valuable tool 
that could be employed in applications requiring a de- 
vice capable of preparing specific quantum mechanical 
states needed for quantum information processing. In 
this paper, we model the full 3D dynamics of this coher- 
ent quantum process, treating the light fields within the 
semiclassical approximation. 

Figure [l] illustrates the generation of a rotating cloud 
of atoms, which travels at a large velocity with respect to 
its stationary source condensate. In the experiment [5], 
this is achieved by applying to the initial BEG suitably- 
tuned counter-propagating Gaussian (G) and Laguerre- 
Gaussian (LG) laser beams. The applied light is linearly 
polarized and therefore has zero-average spin-angular 
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FIG. 1: (Color online) Visualization of the computed conden- 
sate density (rhs) after the Bragg scattered vortex state, de- 
noted by |1), has spatially separated from the initial ground 
state, |0), atoms. The frames in the left are column densi- 
ties integrated along the z-axis of those two distinct angu- 
lar/linear momentum states. The component |1) is traveling 
to the direction of the negative z— axis and the state |0) is at 
rest in the laboratory. 



momentum. Linear and orbital angular momenta car- 
ried by the photons are consequently transferred to the 
BEG in a coherent two-photon Raman process in which a 
condensate atom absorbs one unit, fi, of angular momen- 
tum and one unit, hk^ of linear momentum from a LG 
photon. Correspondingly, in order to satisfy energy and 
momentum conservation laws, the atom emits a stim- 
ulated photon to the G field acquiring a further hk of 
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linear momentum from the recoiling photon. The linear 
momentum of each G photon is set approximately equal 
to the respective LG photon. Alternatively, the process 
described above may be cast in terms of coherent Bragg 
scattering of atoms from a traveling optical diffraction 
grating [6] [7]. In this picture, the combined effect of 
the counter-propagating laser fields on the BEG may be 
accounted for in the electric-dipole approximation by a 
chiral light-shift potential: 

^Bragg(r,t) = |AGp + |ALGp+2A^ALGCOs(g2; + Acjt + 0), 

(1) 

where q = 2k^ Au is the frequency difference between the 
G and LG beams, and (j){x^ y) is the spatial phase of the 
LG beam corresponding to a 27r phase winding around 
the axis. The radial dependence of the laser fields in 
the x — y plane (r^ = + ?/^) in terms of the intensities, 
Aq and Ai, and the beam waists, ctq and (Jlg, is 

= V^e-^'/^G and Alq = V^^e-^'/^'c. (2) 
The resulting mean-field Hamiltonian operator 

H = -^y^+g\^ir,t)f+ VtrapW + VBragg(r,t) (3) 

has the Gross-Pitaevskii form containing the time- 
dependent potential, l^ragg^ which describes the coher- 
ent atom-photon coupling, and an external harmonic trap 
potential 

Vt.ap(r) = f Ka:2+c.y+c.,V). (4) 

The main computational difficulty in the full numerical 
treatment of a typical experiment of Bragg scattering 
of a BEG is caused by the high-frequency oscillation of 
the optical diffraction potential, requiring an extremely 
fine spatial grid and a small time step in the temporal 
propagation of the wavefunction. The formal solution for 
this Hamiltonian system for short time intervals, St, 

^(r, t^5t) = e-^^(^)^^/^^(r, t) (5) 

involves an exponentiation of a sum of non-commuting 
operators. To tackle these issues, we have developed an 
efficient, scalable parallel code for the solution of the 
time-dependent non-linear Schrodinger equation. The 
method of choice combines a finite-element discrete vari- 
able spatial representation (FEDVR) with a fourth order 
split-operator technique for the time evolution [8j. The 
power of this procedure derives from the sparse repre- 
sentation of the Hamiltonian matrix in combination with 
an underlying polynomial basis that enables an efficient 
computation of the matrix exponentials from a quadra- 
ture formula. In practice, the representation of the gen- 
eralized non-linear and time-dependent potential is diag- 
onal, while the kinetic energy part in each spatial dimen- 
sion is sparse and may be split into a sum of two block 



diagonal matrices. Finally, the code is parallelized apply- 
ing a 2D cartesian communicator virtual topology within 
the message passing interface (MPI) protocol, resulting 
in marked savings in both memory and computational 
time. 

In accordance with the experiment [5 , our sys- 
tem consists of N = 2 X 10^ interacting ^^Na atoms 
placed in an anisotropic trap with harmonic frequen- 
cies, {uJxjUJy^uJz} = 407r X {2,V^, l}Hz. The coupling 
constant, g = Anfi^a/m, is expressed in terms of the 
5— wave scattering length, a = 2.75nm, and the mass, 
m, of an atom. The laser field intensities are related by 
Ai = 4.3 X lO^^Ao, their waists are, ctq = 175/im and 
c"LG = 85/im, and their temporal shape is modeled by a 
square pulse of varying duration. The frequency differ- 
ence between the G and LG fields is 

Auj = Acoo + Suj (6) 

where Aujq = hq^ /2m = 5000 uOz and duo is the detuning 
from the first order linear Bragg scattering resonance. 
Since the chemical potential, ja/h ^ 38 cOz <^ Acjq, the 
mean- field shift, to the resonance frequency is al- 

most negligible. 

Figure [2] displays some typical numerical results show- 
ing side-view column densities of the system at three dif- 
ferent times. Before the laser fields are applied, the con- 
densate is in the interacting ground state of the harmonic 
trap, |0), and the spatial density of particles is smooth as 
shown in Fig. [2|a). As soon as the chiral optical diffrac- 
tion grating is initiated, part of the original condensate 
scatters into a traveling vortex state, The traveling 
|1) and stationary |0) components interfere in their over- 
lap region, producing a fine helix of fringes, not visible 
in the main frame of Fig. [2|b). The fringe spacing is 
determined by the relative velocity, v = hq/m, between 
the two components. The wavelength of these density 
modulations, magnified in the inset of Fig. |2|b), falls be- 
yond current experimental imaging resolution. Finally, 
Fig. [2|c) displays the system shortly after the two com- 
ponents |1) and |0) have spatially separated. The |0) 
state has become slightly tapered towards the direction 
of the motion of the |1) component due to the interactions 
rendering the radial expansion z— dependent as shown in 
Fig. [2|c). At the end of the Bragg pulse, the trap is also 
switched off, allowing the system to expand ballistically. 
In Fig. [2] the pulse duration is chosen so that the result- 
ing two momentum components have equal populations, 
thus implementing an atomic beam splitter. 

In Fig. |3j we display the fraction of particles in state 
|1) as functions of pulse duration. The experimental 
data for two different Bragg frequencies, 97.5 kHz (•) 
and 100 kHz (o), are shown [9 . Due to a systematic 
centre-of-mass motion of the condensate, these frequen- 
cies are Doppler shifted from the center of the Bragg res- 
onance. Each point in the experimental data correspond 
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FIG. 2: (Color online) Snapshots of the system at three differ- 
ent times, t = (0, 0.6, and 1.2) ms (a-c). Frames show column 
densities, f |(/)(r)|^ dx. Part of the frame (b) has been mag- 
nified in the region {z,y} = {[—24, —22], [10, 12]} /xm shown 
in the inset to resolve the high frequency interference struc- 
ture. Bragg pulse duration, r = 77/is, Aq = 500 huz and 
6w = ooz- 



to a separate experiment and they are obtained by ex- 
plicitly measuring the populations in states |0) and |1) 
after they have spatially separated during a few millisec- 
onds of time- of- flight (TOF). To compare the numeri- 
cal solution with the experimental data, we have plotted 
the expectation value of the 2;— component of the total 
angular momentum per particle, {Lz)/N, as a function 
of the pulse duration shown in Fig. [3] for two different 
pulse parameters Aq = 460 hcOz ^Sw = 50 uJz (blue) and 
Ai^ = 950 hujz^ Sw = 175 coz (red). Since each atom in the 
traveling vortex state carries a single unit of angular mo- 
mentum, the computed curves may be viewed as the frac- 
tion of the transferred population to the vortex state |1) 
as a function of pulse duration. We have also computed 
the expectation value of the total linear momentum per 
particle, {Pz)/N^ scaled by q. Such curves are indistin- 
guishable from the ones shown, leading to the conclusion 
that the transfer of linear and angular momenta is fully 
correlated, verifying the quantized transfer mechanism 
between the light and matter fields. During its exposure 
to the light field, the condensate executes Rabi cycling 
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FIG. 3: (Color online) Fraction of particles in the vortex 
state |1) as a function of the pulse duration. Experimental 
data sets obtained for Bragg frequencies 97.5 kHz (•) and 100 
kHz (o) are obtained by separately measuring particle num- 
ber in both states |0) and |1) after TOF 9 . Computed curves 
show the expectation values, {Lz)/N, for different detunings, 
Slj = 50 (jJz (blue) and Su = 175 Uz (red), and the respective 
intensities Aq = 460 hooz and Aq — 950 hooz ■ 



between the ground state |0) and the vortex state |1) due 
to the photon-mediated coupling between them, as is ev- 
ident in Fig. (3] The maximum atom transfer is limited 
by the overlap between the ground state condensate and 
the spatial shape of the optical Bragg grating. The dis- 
crepancy between the experimental data and computed 
transfer curves may be accounted for by experimental un- 
certainties, such as the residual, in-trap, center-of-mass 
motion of the condensate. 

Since the trap is fully anisotropic, the condensate 
ground state is not an angular momentum eigenstate but 
is instead described in terms of a superposition of differ- 
ent angular momentum basis states. In order to ana- 
lyze the content and transfer of angular momentum to 
the condensate, we express its wavefunction, ^(^,^) = 
Xlp m ^Pj^ ex.-p{ipz/h + imO)^ in terms of the different m 
states. Since the momentum spread around the station- 
ary and traveling momentum components is narrow, we 
only consider two discrete linear momentum values p = 
and p = hq. In Fig.|4j we plot the normalized probability 
densities |cp,rnP of the projection of the wavefunction to 
different angular momentum states as functions of pulse 
duration. In the top frame p = 0, and only even m val- 
ues are occupied. Similarly, in the lower frame p = hq, 
and only odd m states are populated. The moving and 
stationary components have odd and even parities in the 
X — y-plane, respectively. As is seen in Fig. [4]^ a), the 
initial ground state already contains a 3% admixture of 
m = ±2 states due to the trap asymmetry. In the course 
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FIG. 4: (Color online) Populations of different angular mo- 
mentum states as functions of pulse duration. Upper panel 
contains the even parity, p = 0, states and the lower panel is 
for odd parity, p = states. Different m values are indicated 
in the picture. The subfigure in top frame shows evolution of 
populations of all m states in a single plot as function of pulse 
duration. The pulse parameters are as in Fig. [2] 

of time the odd and even m states become populated 
symmetrically around the respective m = 1 and m = 
states. Coupling with the laser fields induces the domi- 
nant dipole transitions between m = and m = 1 states 
and the undulations in the side band m states arise due 
to the mean- field interactions, as is verified by a compar- 
ative simulation in which the interactions are neglected. 
In the inset of Fig. [4] we have plotted populations of all 
{p, m} states as functions of time showing the population 
spreading to larger m values. 

In conclusion, by modeling the recent experiment, we 
have numerically demonstrated coherent transfer of or- 
bital angular momentum from Laguerre- Gaussian pho- 
tons to Bose-Einstein condensed atoms in agreement with 
the experimental observations. We have employed a pow- 
erful numerical method, which enables us to calculate the 
full 3D dynamics of this system. Altering the squeezing of 
the trap frequencies permits the manipulation of the rel- 
ative populations of different angular momentum states. 
This method may be used to implement a beam split- 
ter into two topological states displaying the potential 
for generation of robust phase qubits. Generalizing, the 
LG states of light form a complete orthonormal basis set 
which may be deployed to gain an access to an infinite- 
dimensional Hilbert space of quantum states for quantum 
information processing. The underlying principle allows 
for the creation of exotic optical lattice configurations 
and the study of persistent currents in toroidal trap con- 
figurations [To]. The latter is a promising starting point 
to investigate the feasibility of creating macroscopically 
entangled superflow states of opposite circulation. 
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